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We report the constraints on the dark energy equation-of-state w{z) using the latest 'Constitution' 
SNe sample combined with the WMAP5 and SDSS data. Assuming a flat universe, and utilizing 
the localized principal component analysis and the model selection criteria, we find that the ACDM 
model is generally consistent with the current data, yet there exists weak hint of the possible 
dynamics of dark energy. In particular, a model predicting w{z) < —1 at 2: £ [0.25, 0.5) and w{z) > 
— 1 at z £ [0.5,0.75), which means that w[z) crosses —1 in the range of z £ [0.25,0.75), is mildly 
favored at 95% confidence level. Given the best fit model for current data as a fiducial model, we 
make future forecast from the joint data sets of JDEM, Planck and LSST, and we find that the 
future surveys can reduce the error bars on the w bins by roughly a factor of 10 for a 5-w-bin model. 



I. INTRODUCTION 

The apparent acceleration of the universe discovered 
at the end of last century remains an enigma [1, 2], and 
much effort has been made to find consistent explanation. 
This includes the modification or even the 'reinvention' 
of the Einstein gravity on large scales [3-11], and also the 
addition of the exotic energy budget - Dark Energy (DE) 
term to the right hand side of the Einstein equation. The 
Equation- of- State (EoS) w of dark energy, defined as the 
ratio of pressure and energy density, is usually used to 
classify different DE models. For example, the w of the 
simplest DE candidate, the vacuum energy, is a constant 
of —1, while w is generally considered as a function of 
redshift z for models predicting dark energy dynamics 
such as quintessence [12], phantom [13], quintom [14] and 
so forth. 

It is true that we are still ignorant about the nature of 
dark energy, yet the accumulating high precision obser- 
vational data of Supernova Type la (SNe), Cosmic Mi- 
crowave Background (CMB) and Large Scale Structure 
(LSS) make it possible to study dark energy phenomeno- 
logically, i.e., confronting the numerous dark energy mod- 
els to data and narrowing down the dark energy candi- 
dates by falsifying part of the models, which might be the 
best we can do to approach the truth of dark energy. In 
this data-driven investigation, in most cases one needs to 
parameterize w{z) in the first place, i.e. assume an ad hoc 
functional form of w(z), then fit the related parameters 
to data, reconstruct w{z) and make statements according 
to the result of the reconstruction [15]. Since the result 
may more or less depend on the assumed form of w{z), 
and to minimize the artifacts, one needs to choose the 
parametrization with care - it should be physically moti- 



vated, and statistically sound, i.e., using the least number 
of free parameters to obtain the maximum generality. 

The first example of the viable parametrization is the 
widely used CPL parametrization [16, 17], 



w{a) =WQ+Wa- 



(1) 



where wq and Wa are free parameters. It is widely used 
since it has a simple form, clear interpretation - wq is the 
EoS today and Wa denotes the derivative with respective 
to the scale factor a, thus a dark energy dynamics indica- 
tor - and a small number of free parameters. However, the 
simplicity of the form inhabits the CPL parametrization 
to describe the models whose EoS deviates significantly 
from the linear function of a, for example, a w{z) with 
oscillations [18], or with other local features [19]. 

One alternative is to approximate w{z) using the piece- 
wise constant bins [20] , which is much more general than 
the functional parametrizations. This generality allows 
for the high-resolution temporal reconstruction of w{z), 
based on which one can make further model-independent 
studies using a Principal Component Analysis (PCA) 
method [20-27]. For example, it is possible to know how 
many EoS parameters can be well constrained by cur- 
rent/future data regardless of the forms of parametriza- 
tion; where is (are) the sweet spot(s) (the redshift where 
the error on w{z) gets minimized); and we can even con- 
struct the uncorrelated bands to probe dark energy dy- 
namics more explicitly from the cosmological observa- 
tions [21, 25, 28, 29] ^ . Given the constraints on w{z) from 
the current data, it is useful to know to what extent the 
future surveys will tighten the constraints. In this paper, 
we will focus on the Uncorrelated Band-power Estimates 



^ Aside from PCA, other less model-dependent methods have been 
'Electronic address: Gong-bo.Zhao@port.ac.uk developed, such as smoothing the SN data to derive w{z) [30], 

tElectronic address: xmzhang@ihep.ac.cn constraining the binned dark energy density [31], etc. 
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(UBE) of w{z) from the current and the simulated future 
data. 

In the next section, we wih describe the method and 
data we use in detail and in section III we will present our 
main results, and then we will finish with the summary 
and discussion. 



II. METHODOLOGY AND DATA 



A. Constraining dark energy from current 
observations 



To fit to data, we parametrize our universe as: 



(2) 



where cob = flbh^ and uJc = 0,ch? are the physical baryon 
and cold dark matter densities relative to the critical den- 
sity respectively, Qg stands for the ratio (multiplied by 
100) of the sound horizon to the angular diameter dis- 
tance at decoupling, r denotes the optical depth to re- 
ionization, and , As are the primordial power spectrum 
index and amplitude, respectively. We assume a flat uni- 
verse throughout. 

The dark energy EoS parameters are denoted by A", 
and we consider the following two kinds of parametriza- 
tions in this work, 



w{z) 



Sum of the tanh bins, Xi = {wi}; 

Wo+Wa- z/{l + z), Xii ={wo,Wa}; 



(3) 

where Xi and Xu are the collections of the bin parameters 
for dark energy equation of state (see explanations below) 
and the CPL parameters respectively, and the constraints 
on these parameters allow us to reconstruct the evolution 
history of w{z), which might encode the dark energy dy- 
namics. 

For dark energy parametrization Aj, we approximate 
w{z) using the sum of N piecewise constant bins Wi 
localizing in redshift and vary them to fit data. The 
simplest realization of this binning is to use step func- 
tions [20, 25, 28, 29]. However, the resulting discontinuity 
in 'w{z) makes it difficult to handle dark energy pertur- 
bations, which depends on the time derivative of w{z) 
and plays a crucial role in the parameter estimation, and 
should be treated in a consistent way [32-34]. To solve 
the problem of discontinuity, one can use the smooth and 
differentiable functions such as the cubic spline functions 
[21, 22], and the hyperbolic tangent functions [24] for the 
binning. Here we follow the latter and parameterize w{z) 
as, 



{Wi+i - Wi) 



l+tanh( ^ (4) 



where Wi denotes the value of EoS in the ith bin, and 
Zi, stand for the endpoints of the ith bin. Also note 



that ^ is the transition width of two neighboring bins, 
and is set to 5% of the bin width. We have numerically 
checked that the final result is largely independent of ^ 
as long as it describes a sharp, but numerically stable 
transition. For a total of iV-bin parametrization, we ar- 
range the first — 1 bins to be evenly spaced at low 
redshifts (z < 1), and use one wide bin to model w{z) 
at z > 1 since dark energy becomes less and less im- 
portant as redshift increases, and it has been found that 
there is no hope to resolve dark energy dynamics, if any, 
beyond redshift one even using the future Joint Dark En- 
ergy Mission (JDEM) survey [20, 35]. In our numerical 
analysis, we let all the TV bins (including the high-z bin) 
float in all cases, thus our constraint on w{z) at low-z is 
more conservative than that with fixed high-z bin. Note 
that in our notation, = 1 corresponds to the wCDM 
model where w is a constant regardless of redshift, and 
TV = stands for ACDM model. 

It is true that the larger N we use, the higher tempo- 
ral resolution we can obtain. However, the result will be 
severely diluted for large TV. This is not only due to the 
weakness of current data, but also to the huge degeneracy 
introduced. Therefore we need a criterion to determine 
the optimal number of bins, so that we can detect the 
main features from data by using the minimal number 
of bins. This is an issue of Occam's razor, i.e., we don't 
want to introduce unnecessary parameters. To optimize 
the Goodness of Fit (GoF), we search for an optimal 
in the range N G [0, 10] based on the model selection 
criteria which we will describe in detail in Sec. HI. Then 
w{z) can be reconstructed from data. 

However, the interpretation of the reconstructed w{z) 
might be obscured by the correlations among the dark 
energy bins. To eliminate this ambiguity, one seeks for a 
linear transformation W to rotate the original parame- 
ter vector X into a new parameter vector Q = WX so 
that the resultant new parameters, denoted by the g's, 
are physically and statistically meaningful, namely, they 
directly relate to the w bins and have uncorrelated errors. 

This is an eigenvalue/ vector problem, and can be 
solved by applying the linear manipulations on the co- 
variance matrix C of the w bins (after marginalizing over 
the other cosmological parameters): 



{XiXi^) 



(5) 



We then diagonalize the Fisher matrix F = C^^ so that 
F = O^AO = W^W. Here O is the orthogonal ma- 
trix, A is diagonal and W = F^/^. Then the resulting 
transformation matrix W, obtained by absorbing A^/^ 
into O, can make the new parameters g's uncorrelated 
and the weights (rows of W) are nearly positive definite 
and localized in redshift, which means that qi has a one- 
to-one correspondence to Wi (and close to Wi), making 
the interpretation more transparent [21, 25, 28, 29]. To 
make q{z) = —1 stand for ACDM, we re-scale W so that 

its rows sum up to unity, i.e. W^j = F^/"^ /{Y^k ^ik^)- 
Note, however, that the normalization is arbitrary and 
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the physical inferences do not depend on this normal- 
ization as both the parameters g's and their values for 
a particular theoretical model change consistently [21]. 
Before normalization, the variances of the g's are unity 
since 

Cq = (QQ^) = W{XiXi'^)W^ = Fi/2cFi/2 = I (6) 

While after normalization, the covariance matrix be- 

comes (Ag.Ag,) = <5,,/(Ea ^^Y' This pre- 
scription is called the Uncorrelated Band-power Es- 
timate, or Localized Principal Component Analysis 
(LPGA). 

Parametrization Xu is a widely used functional form 
for dark energy EoS, but it needs to be tested whether 
this linear function of scale factor is general enough to de- 
scribe dark energy for the current data. This can be done 
by reconstructing w{z) from parametrizations Xi and Xu 
respectively, and then making a direct comparison. 

Given the set of cosmological parameters P in Eq (2), 
we calculate the obscrvablcs including the luminosity dis- 
tance. GMB and matter power spectra, and the cosmic 
age using our modified version of CAMB^, which is able to 
calculate all the necessary observables for an arbitrary 
w{z). We pay particular attention to the Dark Energy 
Perturbations (DEP), especially when w(z) crosses — 1. 
Our calculation for DEP for general dark energy models 
is based on the work presented in Ref. [33]. We then fit 
to SNe, GMB and LSS observations using a modified ver- 
sion of the Markov Ghain Monte Garlo (MGMG) package 
CosmoMC'^[36] based on the Bayesian statistics. 

The supernova data we use are the recently released 
"Gonstitution" SALT sample [37], and we marginalize 
over the nuisance parameter, which is basically the cal- 
ibration uncertainty in measuring the supernova intrin- 
sic magnitude, in the likelihood calculation. For GMB, 
we use the WMAP five-year data including the temper- 
ature and polarization power spectra [38], and calculate 
the likelihood using the routine supplied by the WMAP 
team^. For the LSS information, we use the Sloan Digital 
Sky Survey (SDSS) Luminous Red Galaxy (LRG) sample 
[39], and marginalize over the bias parameter. Further- 
more, we impose the 1 — a Gaussian priors on the Hub- 
ble parameter and baryon density of ft, = 0.72 ± 0.08 and 
ilijh'^ = 0.022 ± 0.002 from the measurements of Hubble 
Space Telescope (HST) [40] and Big Bang Nucleosynthe- 
sis [41] respectively, and a tophat prior on the cosmic age 
of 10 Gyr < to < 20 Gyr. The total likeHhood is taken 
to be the products of the separate likelihoods C of each 
dataset we used, thus the total is the sum of separate 
from individual observations plus that from the priors 
if we define = ^ 2 \ogC. 



^ http://camb.info/ 

^ http:/ /cosmologist.info/cosmomc/ 

* http://lambda.gsfc.nasa.gov/ 



B. Future Forecast 

Given the current constraints on dark energy, it is use- 
ful to know quantitively how future surveys can improve 
the constraints. Therefore we choose the best fit w{z) 
from current data as a fiducial model, and make a fore- 
cast from the surveys of JDEM, Planck [42] and Large 
Synoptic Survey Telescope (LSST) [43] by employing a 
standard Fisher Matrix technique following [44]. 



1. Observables and Fisher Matrices 

Besides the luminosity distance and the GMB power 
spectra simulated for JDEM and Planck respectively, we 
include the tomographic observables including the spec- 
tra of the Galaxy Counts(GG) and Weak Lensing(Wh) 
from LSST, and all the possible cross-correlations with 
GMB. Mathematically, all the tomographic observables 
we use can be summarized as 

Cr=4n I ^A^/f(fc)//'(fc), (7) 

where is the primordial curvature power spectrum 
and /^'^ (fc) denotes the angular transfer functions. Here 
X,Y e [T, E, G,, £■,], where T, E, d and ej illustrate the 
GMB temperature, E'-mode polarization, the ith redshift 
bin for galaxy counts and the jth redshift bin for weak 
lensing shear respectively. In other words, we consider 
all the possible cross-correlations among GMB, GG and 
WL. 

Given the specifications of the proposed future surveys, 
the tool of Fisher matrix [45] enables us to quickly esti- 
mate the errors on the cosmological parameters around 
the fiducial values. For zero-mean Gaussian-distributed 
observables, such as Cf^, the Fisher matrix is given by 

where Pa{i3) is the a{(3)th. cosmological parameter and 
is the "observed" covariance matrix with elements 
that include contributions from noise: 

Cr^Cr+Nr- (9) 

The expression (8) assumes that all fields ^(n) are mea- 
sured over contiguous regions covering a fraction /sky of 
the sky. The value of the lowest multipole can be esti- 
mated from imin ~ [■'''/(2/sky)], whcrc the square brackets 
denote the rounded integer. 

In general, the noise matrix Nf^ receives the con- 
tributions from both the statistical and the systematic 
errors. For the statistical error, we assume the uncor- 
related Poisson noise on the galaxy overdensity in each 
galaxy bin (Gi) and shear fields (e^), the noise is given 
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FIG. 1: The correlated binned estimates of the EoS from current observational data. The black and red error bars show la 
and 2a uncertainties respectively. The blue dashed line shows ACDM model. 
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FIG. 2: The improved and reduced w.r.t. the ACDM 
model as a function of the number of w bins. The big red stars 
illustrate the optimal model. 
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FIG. 3: The correlation matrices for the w bins for different 
binnings. From upper left to lower right, we show the corre- 
lation matrices for A'^ = 2 to = 10 binnings (high-z bin 
inclusive in all cases, as stated in the main text). 



by [46] 



- 


c "Trms 

= '^u 




= Sij — 




= 0, 



where 7rms is the expected root mean square shear of the 
galaxies, and rij is the number of galaxies per steradian 
in the jth redshift bin. 

Besides the shot noise, we follow [47] and consider 
three classes of the systematics: redshift errors, addi- 



tive errors and multiplicative errors ^. The redshift er- 
rors may stem from three sources: the distortion of the 
total galaxy distribution, z-bias and z-scatter. In our cal- 
culation, we marginalized over 30 Chebyshev expansion 
coefficients describing the distortion to the overall shape 
of the galaxy distribution; one redshift bias parameter 
and one redshift scatter parameter for each redshift bin 



^ Systematics are notoriously difficult to model and predict for the 
forecast, and it can also stem from some physical processes such 
as the baryonic uncertainty studied in [48]. But wc didn't include 
this effect for simplity. 
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to account for the uncertainty of the redshift measure- 
ment of the bins. The additive error can be generated by 
the anisotropy of the Point Spread Function (PSF) and 
they generally present for both galaxy counts and lensing 
shear bins. Following [47], we parametrize the additive 
error as 

{crh^6^ypAfAj{^y (11) 

and choose p = 1, 77 = 0. The fiducial values of the A's are 
chosen to be conservative, (^f )^ = 10~®, [A])"^ = 10~^. 
The multiplicative errors in measuring the shear can be 
parameterized as, 

(57k = (c7).,[i + /. + /,]. (12) 



2. Experiments and C'osmological Parameters 

As mentioned before, the data considered in our fore- 
cast include the SNe observations, CMB power spectra of 
temperature and polarization (T and E), WL, GC, and 
their cross-correlations. We assume CMB T and E data 
from the Planck satellite [42], the galaxy catalogues and 
WL data by the LSST [43], complemented by a futuristic 
SNe data set provided by a future JDEM space mission. 
We strictly follow [44] to set up the the survey param- 
eters. In our forecasts, we use the best fit values of the 
cosmological parameters obtained from current data as a 
fiducial model, and impose a Gaussian prior on the value 
of h from the Hubble Space Telescope (HST) [40]. We as- 
sign a constant bias parameter for each galaxy bin, and 
then marginalize over. 

III. RESULTS 

Before viewing the results for parametrization Aj, one 
might be able to make an 'intuitive guess' based on the 
following reasoning. As elaborated in Sec. II, we attempt 
to 'see' the possible dark energy dynamics by placing 
numerous w bins around —1 to fit data. If the real dark 
energy EoS were —1, then the number of bins (and even 
the method of binning) shouldn't affect the fit much, i.e., 
the Ax^ would have little dependence on the number of w 
bins^. However, if dark energy were dynamical, e.g., there 
exist some local features in w(z) on some scale, the w bins 
wander around —1 trying to find a better fit to data than 
w = —1, and the Ax^ behavior would be strongly depend 
upon the number of w bins. Consider, if the binning is 
coarse, then it might not be able to resolve the features 



In the limit of infinitely many SNe on the hubble diagram, the 
Ax^ is independent of the binning if the underlying physical 
model of dark energy is cosmological constant. 



in w{z) properly, and this might result in the marginal 
improvement on ■ However, if the binning resolution is 
fine enough allowing w{z) to vary on the scale of main 
feature of 'w(z), the bins may gain sufficient freedom to 
capture the principal features in the data set, giving rise 
to sharp drops on x^- It is true that using finer binning 
can in principle improve further, but it is not effective 
to improve the reduced x^j which is an indicator of the 
goodness of fit, and is defined as, 

X?cd = xVi^ (13) 

where v denotes the number of degrees of freedom in 
the fit, i.e., the number of data points subtracted by the 
number of fit parameters defined in Eq. (2). For a reliable 
fit, the Xred should be close to unity, which is the case 
for all our 10 fits shown in Fig 1. For example, for our 
base model (ACDM), the x^ is 3149.06 for 2978 degrees 
of freedom, giving Xrcd ~ 1.057. 

Using a unnecessary large number of bins can hardly 
improve the fit. On the one hand, the super fine bins can 
do nothing more but resolve the unimportant feature in 
the data, or overfit the data by treating the noise as fea- 
tures, so that they cannot improve the x^ drastically. On 
the other hand, the redundant bins will be highly cor- 
related, making it difficult to extract useful information 
from data by performing the global fit even if the MCMC 
algorithm is used. Therefore, there is a trade-off between 
the number of w bins N and the GoF, and it is a necessity 
to seek for the optimal maximizing the GoF. 

Now let's look at the reconstructed w{z) from the cur- 
rent data shown in Fig 1. With respect to the ACDM 
model as a base model, we plot the improvement of the 
X^ in panel (A) of Fig 2 and in Table I. We find that as 
A^ increases, the x^ generally decreases, showing an im- 
provement of the fit with more w bins which is expected. 
Furthermore, the slope of Ax^ varies with A^, namely, 
when N goes from 3 to 5, the x^ drops quickly whereas 
A' < 3 or A^ > 5, the x^ decreases slightly. As A^ in- 
creases from 1 to 3, 'w{z) starts to deviate from —1. But 
this makes the x^ gst only slightly improved, meaning 
that these binnings are too coarse to capture the fea- 
tures properly. When A^ grows to 4, a new pattern ap- 
pears, say, all the w bins show deviations from —1 at 
68% confidence level, and interestingly, the odd bin val- 
ues become greater than —1 whereas the even bins drop 
below the —1 boundary. Compared to the case of A^ = 3, 
the resolution of A^ = 4 allows us to see a new feature, 
say, w{z) < — 1 (2 e [0.33, 0.66]). This new feature might 
explain the 'drop' on the Ax^ plot. As A^ reaches 5, the 
pattern of w{z) resembles that of A^ = 4, namely, 

w{z)>-l, ze5 = [0,0.2] U [0.4, 1.0]; 

w{z) < ~1, zeS. (14) 

but the dynamics is seemingly more pronounced - the 
deviations from —1 of the first three bins are enhanced 
to about 95% confidence level. This makes the Ax^ take 
a nosedive again. When A^ exceeds 5, the w{z) pattern 
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TABLE I: The x^, Xrcd (reduced x^) and AIC, BIG values for various models. For each quantity, the difference with respect to 
the base model ACDM is also listed. 



Model 




Ax' 


X?cd(xl03) 


Ax?cd(xl0^) 


BIG 


ABIC 


AIC 


AAIC 


ACDM 


3149.1 





1057.0 





3197.1 





3161.1 





N = 1 


3147.8 


-1.3 


1056.9 


-0.1 


3203.9 


6.8 


3161.9 


0.8 


N = 2 


3146.8 


-2.3 


1056.9 


-0.1 


3210.8 


13.7 


3162.8 


1.7 


N = 3 


3146.2 


-2.9 


1056.9 


-0.1 


3218.2 


21.1 


3164.2 


3.1 


N = 4 


3143.3 


-5.8 


1056.5 


-0.5 


3223.3 


26.2 


3163.3 


2.2 


N = 5 


3139.7 


-9.4 


1055.6 


-1.4 


3227.7 


30.6 


3161.7 


0.6 


N = 6 


3139.5 


-9.6 


1056.0 


-1.0 


3235.6 


38.4 


3163.5 


2.4 


N = 7 


3139.4 


-9.7 


1056.3 


-0.7 


3243.4 


46.3 


3165.4 


4.3 


N = 8 


3139.3 


-9.8 


1056.6 


-0.4 


3251.3 


54.2 


3167.3 


6.2 


N = 9 


3139.2 


-9.9 


1056.9 


-0.1 


3259.3 


62.2 


3169.2 


8.1 


N = 10 


3139.2 


-9.9 


1057.6 


0.6 


3267.2 


70.1 


3171.2 


10.1 


{wo, Wa} 


3146.7 


-2.4 


1056.9 


-0.1 


3210.7 


13.6 


3162.7 


1.6 



TABLE II: The mean values of the dark energy parameters 
with 68% and 95% C.L. error bars for N = 5 model and for 
the CPL model. For the N = 5 model, the constraints on the 
rotated parameters q's are also listed. 



_n 70+0.12+0.23 

,„„ _2 1+0.49+0.93 
W2 Z.i_o5i_o.79 

W3 omtl'otli 

IIU 41+1.9+2.6 

_1 K+0.88+1.3 
W5 l.O_0.9l_l.4 


„, -n qQ+o.07+a.i3 

1 9+0.14+0.26 
q2 l.^_o.l4-0.24 

q.i U.4D_o.28-0.57 

n. n qo+0.55+0.83 
(74 -U.rf8_o5g_io 

1 1+0.37+0.62 
ys l.l_o.34-0.46 


„,„ _o qn+o. 11+0.23 
wo u.yu_o ji_o 22 


Wa -0.24l°;r,l?:f 



remains, but the gets improved only marginally since 
all the main features have already been identified by 5 
bins and there is little work left for the extra bins to do. 
This is in consistency with our foregoing 'intuitive guess' 
for the case of dark energy with dynamics. 

To quantify the GoF for different N, we need to view 
the reduced x^ versus the number of w bins N shown in 
panel (B) of Fig 2 and in Table I. As we see, the improved 
Xrcd reaches its extremum at TV = 5. This means that us- 
ing five bins is sufficient to find all the important features 
on w{z), and any extra bins arc redundant. They dilute 
the constraints by introducing degeneracies, but can do 
Httle to improve the fit. Therefore = 5 is the optimal 
number of bins we need. 

Note, however, different model selection criteria may 
prefer different models. Here wc consider two other 
widely used alternatives, the Akaike information crite- 
rion (AIC) [49] and the Bayesian Information Criterion 
(BIC) [50]. They both serve as tools to compare different 
models using a likelihood method, yet they base on differ- 
ent statistical arguments. For example, AIC stems from 
the minimization prescription of the Kullback-Lcibler in- 
formation entropy, yet BIC roots in the approximation 
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FIG. 4: Panel (A,B): The reconstructed w{z) and q{z) using 
5 bins from current data. The vertical black and red error 
bars show la and 2a errors respectively. The blue dashed 
line shows ACDM model; Panel (C): The re-scaled window 
function W; Panel (D): 1-D probability distribution of the 
q's. 



of the Baycs factor. The AIC and BIC are defined as, 

AIC = -2 In jC + 2Np (15) 
BIC = -2 In £ + A^p InNu (16) 

where C is the maximum likelihood, and A'p and Nd 
are the numbers of free parameters and that of the data 
points used in the fit, respectively. Viable models arc sup- 
posed to minimize the quantities of AIC or BIC. The AIC 
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and BIC for different models are listed in Table I. which 
reads, 

1. Both AIC and BIC favor the ACDM model; 

2. Model of = 5 is favored by the reduced crite- 
rion, and is mildly disfavored by AIC; 

3. All models except ACDM are strongly disfavored 
by BIC. 

It is not surprising that AIC and BIC resist dynamical 
dark energy models because both AIC and BIC include 
strong penalty terms inhabiting the overfit - reducing 
by introducing redundant free parameters. And if the 
number of data points exceeds ~ 7, which is often 
the case in cosmology, the penalty of BIC for additional 
parameters is stronger than that of the AIC. To be ex- 
plicit, let's suppose one model has AA^p additional free 
parameters compared to ACDM. If these AA^p extra pa- 
rameters can help reduce by 2AA'p, then it is preferred 
to ACDM by AIC. However, to satisfy BIC, the must 
be reduced by three times more, i.e. Ax^ ^ 8ANp. 

However, no matter favored or not by AIC or BIC, 
the dynamical dark energy model for the N = 5 case is 
worth investigating in depth since it reduces the x^ most 
for each degree of freedom, and the reconstructed w{z) 
shows an excellent convergence. But unfortunately, this 
reconstruction is blurred by the correlations among all 
the w bins, making Fig 1 and Eq. (14) hard to interpret. 
The correlation between bins Wi and Wj is defined as. 



Corr(u;i, Wj) 



Cov{wi, Wj) 
cj{wi)a{wj) 



(17) 



where Cov(u;i, Wj) is the covariance between bins Wi and 
Wj, and a{wi) and <y{wj) are standard deviations of Wi 
and Wj respectively. In Fig 3, we plot the correlation 
matrices for the cases of A^ = 2 to A^ = 10. As we can 
see, for the w bins centering at z < 0.5, there exists 
strong correlation between the neighboring bins. The last 
bin {z € [1, 1000]) has little correlation with other bins, 
which explains why the constraint on this bin doesn't get 
diluted even if we have large number of redundant bins, 
as shown in Fig 1. 

To de-correlate the bins, we apply the LPCA procedure 
explained in Sec. II to rotate the w's into g's, and we find. 



> 


-1 (68% C.L.), z 


e [0,0.25); 


< 


-1 (- 95% C.L.), 


z e [0.25,0.5); 


> 


-1 (~ 95% C.L.), 


z e [0.5,0.75); 


> 


-1 (68% C.L.), z 


6 [0.75,1); 




-1.44,-0.73] (68^ 


Vo C.L.), z>l. 



(18) 



So at low redshift z < 0.25 or at high redshift z > 
0.75, q{z) is consistent with the ACDM prediction at the 
95% confidence level. However, in the redshift range [0.25, 
0.75), q{z) crosses the cosmological constant boundary, 



and the error bars on the g's are uncorrclatcd thus free 
of degeneracy by design. 

The detailed result of the LPCA is summarized in Ta- 
ble II and in Fig 4. From Panel (C) in Fig 4, we see that 
the window functions are almost positive and fairly local- 
ized, making the g's directly relate to the w's. This means 
that the q's have almost one-to-one correspondence with 
the original w bins. Therefore we can come to conclu- 
sion that the ACDM model is generally consistent with 
current data, yet Eq. (18) implies some hint, very weak 
though, on the possible dynamics of dark energy. In par- 
ticular, a model predicting w{z) < — 1 if z e [0.25,0.5) 
and w{z) > —1 ii z £ [0.5,0.75), which means that w{z) 
crosses —1 in the range of z S [0.25,0.75), is mildly fa- 
vored. 

Based on the best fit w{z) from current data as the 
fiducial model, we can make the future forecast from 
JDEM, in combination of the Planck and LSST surveys. 
The results are shown in Fig 5. As we can see, the win- 
dow functions are localized as that for the current data, 
and the error bars on w{z) get shrink by roughly a factor 
of 8 and 13, for the cases with systematics and without 
systematics, respectively. Note that including systemat- 
ics changes the correlations among the dark energy bins, 
which gives rise to the corresponding changes in W and 
the g's. Given the constraints on the g's, we can conclude 
that, if the fiducial model we derived from current data 
were true, then even for the most conservative case where 
the full systematics are included, JDEM combined with 
Planck and LSST would be able to detect a roughly lOfi 
deviation of w{z) from —1. 

It is useful to study the dependence of the results on 
the form of dark energy parametrizations. For an illus- 
tration, we compare the result we obtained so far to that 
from parametrization Aii, i.e., the CPL parametrization, 
and the result is summarized in Table I and in Fig 6. As 
shown, we find that the constraints on the CPL param- 
eters from current data are 



Wo - -0.90_o. 11-0.22' - -0.24_o.55-1.2 • 



(19) 



The central values indicate that the 'quintom B' scenario 
is mildly favored, namely, the EoS today w{z)\z=o = 
wq > —1, while EoS in the far past w{z)\z=oo = wg+Wa < 
— 1. This is consistent with the recent published result 
using the 'Constitution' SNe sample [51]. The contour 
plots of Wq and Wa are shown as shaded grey regions in 
the right panel of Fig 6. We can see the best fit model 
from current data lies within the 'quintom B' region, and 
the ACDM model is consistent with current data at la. 
Based on the current best fit model, future surveys has 
the ability to exclude ACDM model at 2(7 and 5(t for the 
cases with and without the systematics respectively. In 
the left panels, we make a direct comparison of the re- 
constructed w{z) from two different parametrizations - 
binned w and the CPL, shown in error bars and shaded 
regions respectively. Interestingly, we find that the N = 2 
result agrees well with CPL (and the x^ for these two 
fits are very similar, see Table I), whereas for the case 
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FIG. 5: Future forecasts from JDEM+Planck+LSST on w{z),q{z) and the window functions. Left panels: Forecasts including 
statistical errors only; Right: Statistical plus systematic errors. The legends are the same as that for panels (A-C) in Fig 4. 




Quintessence 




Quintijm 




redshift z 



FIG. 6: Left: Comparison of the reconstructed 'w{z) using two different parametrizations. The 68%(black) and 95%(red) C.L. 
error bars show the reconstructed w{z) using parametrization Xi; The shaded regions show the w{z) reconstruction using 
parametrization Xu. The inner and outer shades illustrate the 68% and 95% C.L. errors respectively, and central white line 
show the best fit model; Right: Contour plots for wo, Wa- Grey: current data; Yellow: forecast from JDEM+Planck+LSST with 
systematics; Blue: forecast without systematics. For the same color, the dark and light shaded regions denote the 68% and 95% 
C.L. contours respectively. The red star stands for the ACDM model, and the white cross illustrates the current best fit model. 



of iV = 5, there is an apparent discrepancy. This is as 
expected - the results converge if the numbers of free 
parameters in the fit are the same, however, the CPL 
parametrization does not have enough freedom to resolve 
the local details of w{z), resulting in the failure to cap- 
ture the 'dip' and 'bump' happening at 0.25 ^ z < 0.75. 



IV. DISCUSSION AND SUMMARY 

In this work, we have investigated the constraints on 
the general form of the equation-of-state of dark energy 
from the latest SNc, CMB and LSS data. We utihze a 
model-independent strategy - redshift binning plus PCA 



- to extract information from data as much as possible. 
Starting from the most general parametrization - the w 
binning using the smooth tanh bins, we constrain the w 
bins using a MCMC algorithm while paying particular at- 
tention to the consistent implementation of the dark en- 
ergy perturbations. We repeat this procedure for different 
binning scheme and investigate the goodness of fits us- 
ing three different model selecting criteria - the reduced 
X^ AIC and BIG. While AIC and BIG strongly favor the 
ACDM model, the 5-w-bin model gives the maximum re- 
duced x^, and we find a convergent evolution trend of 
w{z) when > 4. 

Wc choose the model of A'' = 5 and rotate the w bins 
into the q bins to eradicate the correlations using a LPGA 
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method, and we found that at low redshift z < 0.25 or 
at high redshift z > 0.75, q{z) is consistent with —1 at 
the 95% confidence level. However, in the intermediate 
redshift range [0.25, 0.75), q{z) crosses the cosmological 
constant boundary. Since the q's have almost one-to-one 
correspondence with the original w bins by design, we can 
draw the conclusion that the ACDM model is generally 
consistent with the current data, yet there exists some 
weak hint of the possible dynamics of dark energy. In 
particular, a quintom model predicting w{z) < — 1 if z g 
[0.25,0.5) and w(z) > -1 if z £ [0.5,0.75), which means 
that w{z) crosses —1 in the range of 2 e [0.25,0.75), is 
mildly favored. 

Note that the measurement of the luminosity distance 
is crucial to study dark energy dynamics, thus our re- 
sults are sensitive to the choice of the SN data set and 
analysis. The apparent dynamics of dark energy we found 
using the published 'Constitution' SNe sample might be 
physical, but there is another possibility of the artifacts 
in the SNe data analysis. Therefore we have planned a 
careful study of the effect of SN sample on dark energy 
constraints. 

Given the best fit model from current data, we make 
a forecast from the upcoming/future surveys of JDEM, 
Planck and LSST, and wc find that the future data are 
able to shrink the error bars on the dark energy bins by 
roughly a factor of 10, which is promising to find the 



smoking gun of the dark energy evolution. 

Note added: After completion of this work, we no- 
ticed another analysis on w{z) using different dataset and 
method [22] . Although both our work and Ref. [22] claim 
that the cosmological constant is favored at 95% CL., 
Ref. [22] found less preference of dark energy dynamics. 
This difference may stem from the difference in data se- 
lection and analyzing prescription, and it is worth further 
investigating. 
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